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We study the time evolution of a single spin coupled inhomogeneously to a spin environment. Such 
a system is realized by a single electron spin bound in a semiconductor nanostructure and interacting 
with surrounding nuclear spins. We find striking dependencies on the type of the initial state of 
the nuclear spin system. Simple product states show a profoundly different behavior than randomly 
correlated states whose time evolution provides an illustrative example of quantum parallelism and 
entanglement in a decoherence phenomenon. 



I. INTRODUCTION 

The interest in electron spin dynamics in semiconduc- 

[ tor structures has remarkably increased in the recent 
years generating the emerging field of spintronics^'^. This 
key word summarizes the efforts to use the spins of quan- 
tum objects rather than or in combination with their 
charge for information processing, or, even more ambi- 
tious, for quantum information processing. Meanwhile 
several proposals for quantum information processing us- 

I ing (electron or nuclear) spins have been put forward^"^. 
In quantum information processing the coherence of 

■ quantum bits is crucial. This issue becomes particu- 
larly important in solid state systems where the qubits 
are ususally assumed to be affected by much more and 
much stronger perturbating influences than in other ex- 
perimental setups related to quantum information pro- 

\ cessing in various other fields of physics such as atomic 
physics, quantum optics, or NMR experiments. An im- 
portant advantage however solid state systems have is 
that they offer the possibility scalability once individual 
qubits and elementary gate operations between them are 
established. Such a perspective is usually not given in 
other quantum computation scenarios. 

[ Motivated by these developments we investigate in this 
work the time evolution of single spin S that is cou- 

I pled inhomogeneously to a non-interacting environment 
of other spins. A natural realization of such a system is 

] given by the spin of a single electron in a semiconductor 
quantum dot interacting with surrounding nuclear spins 
via hyperfine coupling®. Alternatively one can think of 
an electron bound to a phosphoros atom implantatcd into 
a silicon matrix^, or of other combinations of impurities 
and host materials such as Si in a GaAs or Ge matrix. 
The time evolution of an electron spin under such circum- 

[ stances is of fundamental interest in its own right and of 
particular relevance to the quantum computing proposal 
of Ref.^'^. In fact, very recently a series of studies of 
electron spin dynamics related to the present one has 
appeared^"^^. Here we build on recent work in Ref.^^, 
where the dynamics of an electron spin due to the hy- 
perfine interaction with nuclear spins in a semiconductor 
quantum dot was investigated. This scource for spin de- 
cay can be assumed to be the dominant one in a quantum 



dot geometry where other mechanisms induced by spin- 
orbit interaction are believed to be suppressed, although 
this issue has not been clarified entirely yet-'^^. 

Our approach here is based on numerical simulations 
of the full quantum mechanical spin dynamics in suf- 
ficiently small systems. These studies complement the 
earlier work in Ref.^'^, where approximation- free analyt- 
ical results where obtained for the case of a fully polar- 
ized system of nuclei. For a more general initial condi- 
tion such as an unpolarized nuclear spin system low-order 
time-dependent perturbation theory was employed which 
unfortunately suffers from divergent terms in higher or- 
der. 

In our numerical simulations we observe a decay of the 
electron spin as measured in terms of the expectation 
values of its components. We study this phenomenon 
as a function of the initial nuclear polarization, and the 
type of the initial state of the nuclear spin system. We 
compare the time evolutions of initial states where the 
nuclear system is in a simple tensor product state with 
the situations where the initial nuclear state is randomly 
correlated. A major results is that the time evolution de- 
pends very significantly on the type of initial state of the 
spin environment. The time evolution of simple tensor 
product states can be quite individual, while randomly 
correlated (and therefore highly entangled) states show 
very reproducible dynamics that mimics the average over 
the time evolutions of all possible tensor product states. 
This observation is an example of quantum parallelism in 
a decoherence phenomenon. 

The spin decay is adccompanied by the generation of 
quantum correlations between the electron spin and the 
nuclear spins, illustrating a general concept of quantum 
information theory where the decoherece of a quantum 
bit (here the decay of the electron spin) is viewed as the 
result of the generation of entanglement (i.e. quantum 
correlations) between the qubit and its environment. We 
quantify this entanglement using well-established meth- 
ods and concepts of quantum information theory. By this 
we also hope that studies of this kind will faciliate fruit- 
ful interactions between the communities of solid state 
physics and quantum infromation. 

Finally we compare the results for the full quantum 
mechanical dynamics with simulations of a classical spin 
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model that arises as the classical limit of the underlying 
Hamiltonian. This comparision shows that the spin de- 
cay observed in the quantum system depends crucially 
on two properties of the system: (i) the inhomogeneity 
of the hyperfine contact interaction induced by the spa- 
tial variation of the electron wave function, and (ii) the 
quantum mechanical nature of the dynamics allowing for 
nontrivial correlations (entanglement) between the elec- 
tron spin and the nuclear spins. 

This paper is organized as follows: In section II we 
briefly describe the details of our modeling and techni- 
cal aspects of our numerical simulations. In section III 
we report on our numerical results. We discuss the role 
of different initial conditions for the nuclear spin system, 
and the connection between decoherence and the gener- 
ation of entanglement observed in our simulations. We 
close with conclusions in section IV. 



II. THE MODEL 

We consider single spin S which is coupled inhomoge- 
neously to a non-interacting environment of other spins, 

i 

The coupling is inhomogcncous since the constants Ai 
vary among the environment spins /j. For an electron 
spin residing in semiconductor quantum dot the coupling 
constants Ai are given by Ai = Auo|^'(r*;)P whore A is 
an overall coupling parameter and Vq the inverse density 
of nuclei in the material, ^{fi) is the electron envelope 
wave function at location f^. This factor induces a spa- 
tial dependence of the cupliong constants Ai which is 
crucial for the spin dynamics. For simplicity we will con- 
sider in the following nuclear spin of length / = i in a 
spherical quantum dot In our simulations a given num- 
ber N of nuclear spins is contained in a sphere of radius 
R = {3N/4Trnoy^^ where no = I/vq is the density of 
nuclei. The electron wave function is given by 

where the parameter a describes the confinement of the 
electron in the dot. In the following we shall use a — 2 
and the material parameters of gallium arsenide with 
no = 45.55nm~^. Therefore a typical quantum dot con- 
tains about TV = 10^ nuclei. For the alternative scenario 
of an electron bound to a ^^P in silicon the number N 
of nuclear spins effectively coupled to the electron spin 
is smaller. The Bohr radius of the hydrogen-like elec- 
tron state is about 3nm, and with the lattice constant of 
silicon and the natural abundance of ^^Si this leads to 
values for N of about a few hundred. 

To mimic their sperical distribution also in systems of 
smaller size used in our simulations we choose the radial 



coordinate of the i-th nucleus according to fj = [3(i — 
l/2)/(47rno)]i/3 with i ranging from 1 to TV. 

The Hamiltonian (1) does not include the direct dipo- 
lar interaction between nuclear spins. This interaction is 
by orders of magnitude weaker than the scale A of the 
hyperfine coupling which is of order 10~^eV in GaAs^. In 
this material the characteristic time T2N for the nuclear 
spin decay due to dipolar interaction is of order 10~^s 
while the time scales considered in this work will be at 
least two orders of magnitude smaller. We also mention a 
recent interesting numerical study by Dobrovitski et al}^ 
on spin dynamics stressing the role of entropy. There a 
central spin is coupled inhomogeneously to an essentially 
non- interacting spin environment, where, differently from 
the present study, an Ising-like coupling was used. To 
allow for nontrivial dynamics the authors of Ref.^* intro- 
duced a magnetic field perpendicular to the ^-direction 
of the Ising coupling. 

The model (1) was specifically studied recently in 
Rcf."'^'', where approximation- free analytical results were 
obtained for the case of a fully polarized system of nu- 
clei. For a more general initial conditions such as an un- 
polarized nuclear spin system low-order time-dependent 
perturbation theory was employed which unfortunately 
suffers from divergent terms in higher order. In this work 
we choose a different route and perform finite-size exact 
diagonalizations from which we obtain the full time evo- 
lution. 

Since our Hamiltonian conserves the total spin J = 
S + it is convenient to work in a subspace of given 
having a dimension of 

( J. ) • (3) 

To obtain the time evolution of the total spin system with 
the initial state lying in a subspace with a given value 

we diagonalize the Hamitonian within this subspace 
and compute the time evolution of certain expectation 
values from the eigensystem data. For initial states hav- 
ing non-zero overlap in several subspaces with different 

the time evolutions obtained in the different spaces 
have to be superimposed. The fact that the full eigen- 
system of the Hamiltonian is required in this procedure 
is different from most other numerical investigations of 
spin systems where, for reason of the physical questions 
being investigated, it is sufficient to concentrate on the 
ground state and some low- lying excitations. In our case 
we need the full eigensystem and are therefore restricted 
to system sizes N and values of yielding dimensions of 
not more than a few thousands. However, as we will ex- 
plain in more detail below, our findings neither depend 
on the restriction to a certain value of nor on the 
specific coupling constants induced by the quantum dot 
geometry. In fact, qualitatively the same results are ob- 
tained when working in the full Hilbert space, or if the 
coupling parameters are chosen at random from a uni- 
form distribution. On the other hand it is essential that 
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the coupling is inhomogeneous leading to a time evolu- 
tion which is for all practical time scales aperiodic and in 
this sense irreversible. The case of homogeneous coupling 
with all Ai being the same is readily solved analytically 
and generates periodic dynamics with recurrence time 
T = AnNn/A. 

III. RESULTS FOR THE TIME EVOLUTION OF 
SPINS 

In the time evolutions to be discussed below the initial 
state will always be a simple direct product of the state 
of the nuclear or environmental spin system and the elec- 
tron spin with the latter pointing downwards along the 
z-direction. Therefore the nuclear spins and the elec- 
tron spin are initially uncorrelated. For the nuclear spin 
system itself we consider two types of initial conditions 
which give rise to significantly different time evolutions. 

A. Product states versus randomly correlated states 

We investigate two types of initial states for the nu- 
clear spin system which differ crucially in their corrrela- 
tion properties and, as a we shall see below, as a conse- 
quence of this also in their time evolution: (i) The nu- 
clear spins are initially in a simple tensor product state. 
If is fixed to a certain value such a state consists 
of eigenstates of If for each nuclear spin i. If this re- 
striction is not applied tensor product states consisting 
of more general spin-coherent states are possible, (ii) 
The nuclear spin state Ix^v) is initially a linear super- 
position, \xn) = J2t'^t\T), where the sum goes over 
all tensor product states |T) consisting of eigenstates /f , 
i G {1, . . . , N}, and is, for fixed J^, restricted to the ap- 
propriate subspace. The coefficients ax in this entangled 
pure state are subject to a normalization condition and 
chosen either at random or coherently (for example, can 
have the same phases). 

As we shall sec shortly, a single tensor product state 
on the one hand and a randomly correlated nuclear state 
on the other hand generate strikingly different time evo- 
lutions for the electron spin. Fig. 1 shows numerical data 
for time-evolved expectation value {S^{t)) for an initially 
randomly correlated system and different degrees of its 
polarization (characterized by J^). In all cases, {S^{t)) 
decreases in magnitude. With decreasing polarization the 
decay becomes more pronounced, and the oscillations ac- 
companying this process get suppressed. Note that it is 
the decay of the envelope in these graphs but not the 
fast oscillations itself that signals the decay of the spin. 
The distance between two neighboring maxima of the 
oscillations can depend slightly on the initial state and 
the coupling constants in the Hamiltonian. However, a 
good estimate for this effective period is usually given by 
T = A-kTi/A since A/2 is an estimate (neglecting quantum 



fluctuations) for the width of the spectrum, i.e. the dif- 
ference between the largest and the smallest eigenvalue 
of the Hamiltonian. 

When the nuclear spin system is initially in a randomly 
correlated state the time evolution of {S^{t)) is very re- 
producible in the sense that it depends only very weakly 
on the particular representation of the initial random 
state. This is illustrated in Fig. 2 where the results of 
different initial random configurations are compared for 
two different system sizes and degrees of polarization. 

This behavior of randomly correlated initial states 
sharply contrasts with the time evolution of simple ten- 
sor product nuclear spin state. The upper two panels of 
Fig. 3 show the time evolution of the electron spin for two 
initial tensor product states. In the lower right panel we 
compare this data with the time evolution of a represen- 
tative of the randomly correlated initial condition. In the 
former case the time evolution depends significantly on 
the concrete initial tensor product state, and the decay 
of the electron spin occurs typically clearly more slowly 
than in the case of an initially randomly correlated nu- 
clear spin system. 

In the lower left panel of Fig. 3 we show the time evolu- 
tion of the electron spin averaged over all nuclear tensor 
product states^'^ Comparing the two lower panels one 
sees that this data is very close to the time evolution of a 
randomly correlated state. This observation is also made 
for other system sizes and degrees of polarization and 
constitutes an example of quantum parallelism^^: The 
time evolution of each initially uncorrelated (and there- 
fore classical-like) nuclear state is present in the evolu- 
tion of a linear superposition of all such states. In other 
words, the time evolutions of all uncorrelated classical- 
like states are perfomed in parallel in the time evolution 
of the randomly correlated state. An experimental con- 
sequence of this observation is that if the electron spin 
dynamics would be detected on an array of independent 
quantum dots one could not distinguish whether the nu- 
clear spin system in each dot was initially randomly cor- 
related or in an uncorrelated tensor product state. In 
other words, the spin dynamics of a randomly correlated 
pure state of the nuclear system in a single dot cannot be 
distinguished from a mixed state of an ensemble of dots. 

The observation that the time evolution of a randomly 
correlated state quite closely mimics the average over all 
tensor product initial conditions relies on the cancellation 
of off-diagonal terms a^ar' {^,T\S{t)\ JJ.,T'), T ^ T', 
due to the randomness in the phases of the coefficients 
ttT- In this sense our system has a self- averaging prop- 
erty. This can be checked explicitly by reducing this 
randomness. The left panel of Fig. 4 shows the time 
evolution of a randomly correlated state where the am- 
plitudes aT are restricted to have a non-negative real 
and imaginary part. This time evolution turns out to be 
similarly reproducible as before, i.e. it does not depend 
on the concrete realization of the initial random state, 
but it is clearly different from the former case since the 
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cancellation of ofF-diagonal contributions is inliibitcd^^ . 
For comparison we show in the right panel data where 
the amplitudes in the initial nuclear spin state have a 
random phase but are restricted to have the same modu- 
lus. Here the proper averaging process takes place again. 
The results described so far were obtained in certain sub- 
spaces of and for the form of coupling constants Ai 
as induced by the quantum dot geometry. However, our 
findings do not depend on these choices. We have also 
performed simulations were the initial state has overlap 
in the full Hilbert space. For a randomly correlated ini- 
tial nuclear spin state the only difference is that now 
also transverse components (5'^(i)), {S^{t)} of the elec- 
tron spin evolve. However, these are tiny in magnitude 
and oscillate around zero. For an initial tensor product 
states these transverse components can become sizable, 
and the time evolution again strongly depends on the 
concrete initial tensor product state. Moreover, as men- 
tioned earlier, the exact form of the coupling constants is 
also not crucial as long as they arc sufficiently inhomo- 
geneous. For instance, we obtain qualitatively the same 
results if we choose the coupling parameters randomly 
from a uniform distribution. 

We also note that coupling a magnetic field to the elec- 
tron spin has only a quantitative influence on our results. 
Here again the time dependence of tensor product initial 
nuclear state is very individual, while a randomly cor- 
related states gives very reproducible results that mimic 
closely the average over tensor product states. 

B. Decoherence and the generation of entanglement 

In circumstances of quantum information processing 
the decay of a qubit is ususally viewed as some 'deco- 
herence' process due to the environment attacking the 
quantum information. As seen above, the spin decay is 
generically slower if the spin environment is initially in 
a uncorrelated state. This finding suggests that it is ad- 
vantageous for protecting quantum information to disen- 
tangle the environment that unavoidably interacts with 
the qubit system. 

A 'decoherence' process of the above kind can be 
viewed as the generation of entanglement between a 
qubit and its environment. The system investigated 
here provides an illustrative example for this statement. 
The entanglement in the total state \^{t)) between the 
central electron spin and its environment can be mea- 
sured by the von-Neumann entropy of the partial den- 
sity matrix where either the electron or the environment 
has been traced out from the pure-state density matrix 
Tracing out the nuclear system we have 



Pel{t) = 




This matrix has eigenvalues \± = 1/2 ± and the 

measure of entanglement reads I ^(t))) = — A+logA+ — 



A_logA_. Thus, the formation of expectation values 
\{S{t))\ + 1/2 (or, in the case of fixed J^ just |(5'^(t))| + 
1/2), is a manifestation of the entanglement between the 
electron spin and the nuclear spin system. The maximum 
entanglement, E = log 2, is achieved if the electron spin 
has decayed completely as measured by the expectation 
values of its components, {S{t)) = 0. The generation of 
quantum entanglement between the electron spin and the 
nuclear spin system signalled by a reduced value of {S{t)) 
is a main and crucial difference between the quantum 
system studied here and its classical 'counterpart' de- 
scribed by a system of Landau-Lifshitz equations. These 
equations can be obtained from the Heisenberg equations 
of motion for the quantum system, dS/dt = i[H,S]/h, 
dli/dt = ilH^Iil/h, by performing expectation values 
of both hand sides within spin-coherent states and as- 
suming that the expectation values of all operator prod- 
ucts factorizes to products of expectation values. This 
procedure becomes exact in the classical limit^''. The 
resulting equations do not contain operators any more 
but just describe the dynamics of three-component vec- 
tors (classical spins) of fixed length. We have performed 
simulations of such a classical spin system by solving 
the Landau-Lifshitz equation via the fourth-order Runge- 
Kutta scheme. As a result, the central classical spin per- 
formes an irregular chaotic motion which does not show 
any similarity to the results for the quantum spin- 5 case. 
In particular all qualitative features of quantum effects 
such as the generation of entanglement (signalled by a 
decay of spins as measured by their expectation values) 
are not present in such a time evolution. Therefore the 
Landau-Lifshitz equation provides only a rather poor de- 
scription of the underlying quantum system. 

Let us now briefly discuss how the different initial con- 
ditions can be prepared experimentally. A tensor prod- 
uct initial state can be produced by applying a magnetic 
field and having the underlying crystal lattice at a tem- 
perature high enough such that spin-lattice relaxation 
processes to the nuclear spins are efficient. These interac- 
tion with the phonon environment will effectively perform 
projection-type measurements on each spin and force the 
system to be in a state close to a tensor product of nuclear 
states pointing in each of the two direction along the field 
axis. Another possibilty is the use of all-optical NMR 
techniques as described in Ref.^^ . A randomly correlated 
nuclear state on the other hand can be achieved by cool- 
ing down the lattice to temperatures where phonon pro- 
cesses are suppressed. Then the highly anisotropic and 
long-ranged dipolar interaction will produce a sufficiently 
'disordered' state with a highly irregular pattern of am- 
plitudes when expressed in the tensor product basis, as 
we have confirmed by explicit simulations of a system of 
8 nuclear spin placed on the edges of a cube. The highly 
correlated (or entangled) character of these states can 
be detected by following the individual nuclear spins in 
terms of their expectation values \{Ii{t))\. This quantity 
decays from its initial value of 1/2 (in a tensor product 
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initial state) on a time scale determined by the dipolar 
interaction^"^ to values typically close to zero. According 
to the entanglement measure E discussed above this in- 
dicates strong entanglement between each nuclear spin li 
and its environment of all other nuclear spins. In both 
cases the initial state of the full system can be prepared 
by injecting the electron to the quantum dot from an 
external lead, or the electron state can be prepared by 
cooling in a magnetic field and ESR techniques. 

We finally consider the nuclear spin correlator C{t) = 
(/^(t)/^(0)), / — ^^li, which can be measured directly 
by local NMR like measurements such as magnetic res- 
onance force microscopy^^. In a subspace of given 
and the electron spin pointing downwards initially this 
quantity reads C{t) = {J'- - (S'^(t)))(J^ + 1/2). A realis- 
tic initial state will have its dominant weight in a series 
of subspaces with neighboring centered around some 
value. Then the time evolution of {S^{t)) is very similar 
in these subspaces, and the dynamics of the total nuclear 
spin can be mapped out by measuring the electron spin 
and vice versa. 



IV. CONCLUSIONS 

In summary we have studied the dynamics of a sin- 
gle spin coupled inhomogeneously to a spin environment. 
As a main result the time evolution depends sensitively 
on the type of initial state of the spin environment. 
While the time evolution of simple tensor product states 
can be quite individual, randomly correlated states show 
very reproducible dynamics that mimics the average over 
the time evolutions of all possible tensor product states. 
This observation constitutes an example of quantum par- 
allelism in a decoherence phenomenon. This effects is 
clearly seen for all finite system sizes studied here and 
can therefore be expected to be also present in realistic 
quantum dot systems containing about N = 10^ nuclei, 
and also in the thermodynamic limit N oo. 

The decay of the single spin in terms of its expectation 
values is due to the formation of entanglement between 
this spin and its environment. Since this decay is gener- 
ally slower if the spin environment is initially in a simple 
tensor product state (i.e. no entanglement among the 
environmental spins) our results suggest that it is advan- 
tagous for protecting quantum information to disentan- 
gle the environment. We expect this result to be of quite 
general nature, i.e. it should also be valid for other sys- 
tems consisting of some central quantum object coupled 
to a bath of other quantum degrees of freedom. 
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time [h/(2jtA)] 

FIG. 1. The time evolution of the electron spin in a system 
of A'" = 14 nuclear spins for different degrees of polarization 
of the randomly correlated nuclear system and coupling con- 
stants induced by the quantum dot geometry. In the top left 
panel the nuclear spins are fully polarized in the initial state 
with the electron spin pointing opposite to them (J^ = 13/2). 
In the following panels the number of flipped rmclear spins in 
the initial state is gradually increased. The case of an initially 
fully unpolaxized (but randomly correlated) nuclear system is 
reached in the bottom right panel {J" = —1/2). Here and in 
the following we take spins to be dimensionless, i.e. measured 
in units of %. 
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FIG. 2. Upper panels: {S'{t)) for a system of size A' = 19 
being initially in a randomly correlated nuclear spin state in 
the subspace with = 7. The two panels represent two 
different randomly chosen initial conditions. Lower panels: 
Analogous data for A = 14 and a completely unpolarizcd nu- 
clear spin system (J^ = —1/2). In both cases the simulation 
data does practically not depend on the initial condition. 
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FIG. 3. Upper panels: time evolution of the electron spin 
{S^{t)) for a system with 14 nuclear spins being initially in an 
uncorrelated tensor product state in the subspace = 9/2. 
The oscillation period and the time scale of the decay are con- 
sistent with the period T = 4iv?i/A and the scale UN/ A iden- 
tified in^''. Lower left panel: data of the same type as above 
but averaged over all possible uncorrelated initial states with 
= 9/2. Here again, the time scale of the decay is consis- 
tent with the scale h\/N /A identified in'^^. Lower right panel: 
{S^{t)) for the same system being initially in a randomly cho- 
sen correlated state. 
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FIG. 4. Time evolution of (S'^(t)) for two types of initially 
randomly correlated nuclear spin states. In the left panel the 
amplitudes olt are restricted to have iiou-iK^gativo real and 
imaginary part, while in the right panel they have all the 
same modulus but completely random phases. 
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